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We consider many-body systems with a global U(l) symmetry on a class of lattices with 
the (fractal) dimensions D < 2 and their zero temperature correlations whose observables 
behave as a vector under the U(l) rotation. For a wide class of the models, we prove that 
if there exists a spectral gap above the ground state, then the correlation functions have 
a stretched exponentially decaying upper bound. This is an extension of the McBryan- 
Spencer method at finite temperatures to zero temperature. The class includes quantum 
spin and electron models on the lattices, and our method also allows finite or infinite 
(quasi) degeneracy of the ground state. The resulting bounds rule out the possibility of 
the corresponding magnetic and electric long-range order. 



Spectral Gap and Decay of Correlations 



2 



1 Introduction 

As is well known, low dimensional systems show large fluctuations for continuous symme- 
try. The most famous result is the Hohenberg-Mermin- Wagner theorem [1] which states 
that the corresponding spontaneous magnetizations are vanishing at finite temperatures 
in one and two dimensions. Since their articles appeared, their method have been applied 
to various systems [2] including classical and quantum magnets, interacting electrons in a 
metal and Bose gas. 1 The theorem was extended to the models on a class of generic lattices 
with the fractal dimensions D < 2 by Cassi [4]. In a stronger sense, it was also proved 
for a class of low- dimensional systems that the equilibrium states are invariant under the 
action of the continuous symmetry group [5, 6]. Even at zero temperature, the same is 
true [7, 8, 9] if the corresponding one- or two-dimensional system satisfies conditions [10] 
such as boundedness of susceptibilities. 2 

About the corresponding long range correlations, Fisher and Jasnow [11] proved clus- 
tering properties of two point functions by using the Bogoliubov inequality. See also ref. [6]. 
McBryan and Spencer [12] obtained a better decay for two point correlations of classical 
spin systems. Their method has been applied various classical and quantum systems. The 
resulting upper bounds for the correlations decay by power, exponential or stretched ex- 
ponential laws [13, 14, 15, 16, 17, 18] and rule out the ordering at finite temperatures in 
(fractal) dimensions D < 2. But a zero temperature analogue of the McBryan- Spencer 
bound has not yet been obtained. 

On the other hand, since Haldane [19] predicted a "massive" phase in low dimensional, 
isotropic quantum systems, many examples have been found to have a spectral gap above 
the ground state and exponentially decaying correlations in the ground state as initiated 
by [20]. See also related articles [21, 22, 23]. These examples raise the following ques- 
tion: Consider a generic low dimensional system with a continuous symmetry and its zero 
temperature correlations whose observables show a non-trivial representation under the 
action of the continuous symmetry group. Then a spectral gap above the ground state 
implies (stretched) exponential decay of the correlations? And if so, the corresponding 
upper bound can be obtained by the McBryan-Spencer method? 

In this paper, we address this problem. We consider the correlation functions whose 
observables behave as a vector under a U(l) rotation. In order to estimate the decay of the 
correlation functions, we extend the McBryan-Spencer method to zero temperature under 
the assumption that there exists a spectral gap above the ground state. As a result, we 
prove that the correlation functions have a stretched exponentially decaying upper bound. 
This method covers a wide class of many-body systems with a global U(l) symmetry on 
a class of lattices with the (fractal) dimensions D < 2. We stress that this method is also 
an extension of the Combes-Thomas method [24] for Schrodinger operators to many-body 
systems in statistical mechanics. In the next section, we will give the precise definition 
of the class of the fractal lattices and describe the results for two typical examples, the 

1 For a mathematically rigorous treatment for the unbounded operators, see ref. [3]. 

2 Since a single spin shows the spontaneous magnetization at zero temperature, the absence of the 
spontaneous symmetry breaking implies that the strong fluctuations due to the interaction destroy the or- 
dering and lead to the finite susceptibilities. In other words, one cannot expect the absence of spontaneous 
symmetry breaking at zero temperature in a generic situation. 
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Heisenberg and the Hubbard models. The proof of the main results will be given in 
Section 3. 

2 Models and results 

We begin with defining the class of the (fractal) lattices which we consider in this paper. 
The class of the lattices is the same as in [18]. See also [25] for fractal lattices and models 
on the lattices. 

Consider first a connected lattice A = (A s , A 6 ), where A s is a set of sites, i, j, k,£, . . ., 
and A?, is a set of bonds, i.e., pairs of sites, {i,j},{k,£},.... If a sequence of sites, 
io,ii,i2, ir, satisfies {i n -i,i n } G A5 for n = 1,2, ...,r, then we say that the path, 
Oo, H,i2, ■ ■ ■ , i r }, has length r and connects io to i r . We define the "sphere", S r (m), 
centered at m e A s with the radius r as 

S r (m) := {£ E A s |dist(£,m) = r}, (2.1) 

where dist(£, m) is the graph-theoretic distance which is defined to be the shortest path 
length that one needs to connect £ to m. Let \A\ denote the number of the elements in the 
set A. We assume that there exists a "(fractal) dimension" D > 1 of the lattice A such 
that the number |5V(m)| of the sites in the sphere satisfies 

sup \S r {m)\ < Cor ' 1 (2.2) 

mgAs 

with some positive constant Cq. 

We consider spin or fermion systems with a global U(l) symmetry on the lattice A 
with the (fractal) dimensions 1 < D < 2. We require the existence of a "uniform gap" 
above the sector of the ground state of the Hamiltonian H A . The precise definition of the 
"uniform gap" is: 

Definition 2.1 We say that there is a uniform gap above the sector of the ground state if 
the spectrum o-(H A ) of the Hamiltonian H\ satisfies the following conditions: The ground 
state of the Hamiltonian H^ is q-fold (quasi) degenerate in the sense that there are q eigen- 
values, £0,1, • • • , Eo,q, in the sector of the ground state at the bottom of the spectrum of H\ 
such that 

AS := max{|E , M - E ^\} — > as |A S | — > 00. (2.3) 

Further the distance between the spectrum, {E 01 , . . . , E q }, of the ground state and the 
rest of the spectrum is larger than a positive constant AE which is independent of the 
volume \ A S \. Namely there is a spectral gap AE above the sector of the ground state. 

Remark: For the special case with q — 1, the ground state is unique. As well known 
examples with q 7^ 1, Majumdar- Ghosh model [21] shows a spectral gap above the degen- 
erate ground state, and the spin-1 antiferromagnetic chain with open boundaries exhibits a 
spectral gap above the fourfold quasidegenerate ground state [22] . In the thermodynamic 
limit |A S | — > 00, we allow infinite degeneracy q = 00 of the ground state as in [26]. 
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Consider first a quantum spin system with a U(l) symmetry on the lattice A. As a 
concrete example, we consider the standard XXZ Heisenberg model on the lattice. The 
Hamiltonian H\ is given by 

H A = H* Y + V A ({Sf}) (2.4) 

with 

Hl Y = 1 E J™(S?S* + S?S]), (2.5) 

{i,j}eA b 

where (Sf, Sf, Sf) is the spin operator at the site i E A s with the spin S — 1/2, 1, 3/2, . . ., 
and jf;J are real coupling constants; V\({Sf}) is a real function of the z-components, 
{Sf} ie \ s , of the spins. For simplicity, we take 

Va({S?})= £ Jf,jS?S* (2.6) 

{i,j}eA b 

with real coupling constants Jf,. We assume that there are positive constants, J^ x and 
^max> which satisfy \Jf^\ < J^ x and \JfJ\ < «/ max for any bond {i, j} £ A&. We stress 
that we can also treat more general interactions in the same way. 

Let P be the projection onto the sector of the ground state. We define the ground-state 
expectation as 

(•••> :=^Tr(...)Po, (2.7) 

where Tr stands for the trace which is over all the spin states. We consider the transverse 
spin-spin correlation, (S+S~) , where Sf := Sf ±iSf. 

Theorem 2.2 Suppose that 1 < D < 2 and that there is a uniform gap AE above the 
sector of the ground state in the sense of Definition 2.1. Then there exists a positive 
constant 7 such that the transverse spin-spin correlation satisfies the bound, 



S^S-) o < Const, exp [- 7 {dist(m, n)} 1 '^ 2 ] , (2.8) 
in the thermodynamic limit \A S \ — > 00. 

Remark: 1. It is easy to extend the result to more complicated correlations such as the 



multispin correlation • • • S+ o S ni • • • S nj / 
2. Recently, the exponential clustering of the correlations was proved for quantum many- 
body lattice systems by Hastings [27, 28] under the gap assumption. This is a non- 
relativistic version of Fredenhagen's theorem [29] of relativistic quantum field theory. 
Combining this exponential clustering with the present result, a better, exponentially 
decaying bound for the ground-state correlations can be obtained [30] for a class of models 
on lattices with the (fractal) dimensions D < 2 and with a certain self-similarity. The 
class includes the translationally invariant regular lattices such as Z. 

For a finite volume |A S | < 00, the Hamiltonian H A of (2.4) commutes with S^ : = 
S^eA s Sf. Using this symmetry, we can block diagonalize the Hamiltonian H A , and denote 
by H\ m the restriction of H\ to the eigenspace 7Ya,m of S^ with the eigenvalue M. Let 
Pq,m denote the projection onto the ground state of Hx,m in the subspace Ha,m- Here 
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the ground state may be (quasi) degenerate in the sense of Definition 2.1. We define the 
ground-state expectation as 



(...) 



0,M 



)P ( 



0,M, 



(2.9) 



where qM is the degeneracy of the sector of the ground state of H\ >M . 

Theorem 2.3 Suppose that 1 < D < 2 and that there is a uniform gap AE above the 
sector of the ground state of H^,m in the spectrum of H^,m with the eigenvalue M of S^ in 
the sense of Definition 2.1. Then there exists a positive constant 7 such that the transverse 
spin-spin correlation satisfies the bound, 



q+ q- 

m n '0,M 



< Const, exp — 7{dist(m, n)} 



l-D/2 



(2.10) 



00. 



in the thermodynamic limit \A S 

As an example of a lattice fermion system, we consider the following Hamiltonian 
[16, 17] on the lattice A with the (fractal) dimension 1 < D < 2: 

H a = ~ E E (kA, c ^ + ttA^) + ^ (K/J) + E Bi • S*, (2.11) 

{i,j}eA b /x=T,i ieAs 

where c| , Cj )At are, respectively, the electron creation and annihilation operators with the 
z component of the spin jj, =t, I, = c\ c ijlx is the corresponding number operator, 

and Sj = (Sf,Sf,S?) are the spin operator given by Sf = ^2^=^,1 c i,/j, (r ^,u c i,f with the 
Pauli spin matrix (a^ u ) for a = x,y,z; G C are the hopping amplitude, V({ni^}) is 
a real function of the number operators, and Bj = (Bf,Bj,B*) e R 3 are local magnetic 
fields. We assume that the interaction V({n itil }) is of finite range in the sense of the graph 
theoretic distance. 

Clearly the Hamiltonian of (2.11) commutes with the total number operator A/a = 



for a finite volume |A S | < 00. We denote by Ha,n the restriction of Ha 



onto the eigenspace of A/a with the eigenvalue A. Let P ,v be the projection onto the 
sector of the ground state of H\ N , and we denote the ground-state expectation by 



(' • V = ^Tr (■ • -)Pa 
qN 



N, 



(2.12) 



where q^ is the degeneracy of the ground state. Assume the existence of a uniform gap 
above the ground state. Then we have 



C m,fM C n,tM/ QN 



< Const, exp — 7{dist(m, n)} 



l-D/2 



and 



C m,'\ C ln,l C ri,'\ c n,l 



0,N 



< Const, exp — 27{dist(m, n)} 



l-D/2 



(2.13) 



(2.14) 



with some constant 7 in the thermodynamic limit |A S 
has the form Bj = (0, 0, Bi), then we further have 



00. If the local magnetic field 



c+ c- 



0,N 



< Const, exp [-7'{dist(m,n)} 1_z?/2 



(2.15) 



with some constant 7'. We remark that, in the latter situation, we can also restrict the 
Hamiltonian Ha,n and the expectation to the subspace with a fixed total magnetization. 
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3 Proof of the main theorems 

We will give a proof only for Theorem 2.3 because the rest of the bounds can be proved 
in the same way. The application to other systems is also straightforward. 

Before proceeding to the proof, let us sketch the idea of the proof and the key tools. 
In the previous work [16], the global quantum mechanical U(l) symmetry was used for es- 
timating the correlation functions of the Hubbard model at finite temperatures. Roughly 
speaking, the strength of the long-range ordering can be measured by twisting the U(l) 
phase locally in the pure imaginary direction. The difference between zero temperature 
and finite temperatures is in their density matrices. Since the projection Po,m onto the 
sector of the ground state can be written in the contour integral of the resolvent, we apply 
the method of [16] to the resolvent (z — Pa)" 1 instead of the Boltzmann weight for finite 
temperatures. It is well known that, for a single-particle resolvent with a classically forbid- 
den energy, the Combes-Thomas method [24] yields the WKB-type tunnelling estimate, 
i.e., the exponentially decaying bound. In the method, twisting the quantum mechani- 
cal phase locally in the pure imaginary direction plays an important role, too, and so the 
Combes-Thomas method is essentially equivalent to the McBryan- Spencer method. In the 
present paper, we use the improved version [31] of the Combes-Thomas method. In order 
to extend the Combes-Thomas method to many-body systems in statistical mechanics, 
we further employ the technique [32, 33] which was developed for many-body systems to 
treat quantities of order of the volume. 

In order to estimate the transverse spin-spin correlation, {S^S~) QM , we introduce 
"gauge transformation" , 

G{oc) := ]J exp^S?], (3.1) 

ieA s 

where a is a real parameter to be determined. For the choice of the function 0j, we 
follow Picco's idea [14], but we modify it in order to obtain a better decay bound for the 
correlation. Let k be a positive parameter satisfying 

D 3D . . 

1 <k< , 3.2 

2 2 2' K ' 



and write R = dist (n, m). We choose the real function 0, on the lattice A s as 

e t = r 1 -°/ 2 x < 



R~ K - 1, for i = m; 

[dist(£, m)/RY - 1, for 1 < dist(£, m) < R; (3.3) 
0, for dist (£,m) > R. 



Note that 

G{a)- 1 Sf G{a) = e ±a6i Sf for i e A s . (3.4) 
Using the relation (3.4), we have 

Tr S+S-P , M = Tr G(«)- 1 5+G(a)G(a)- 1 5+G(a)G(a)- 1 P ,MG(a) 

= e a <*»-«»> Tr S+S-Po^ia), (3.5) 



where we have written as 



P ,m(«) := G(a)- 1 P ,MG'(a). (3.6) 
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This operator has the property, 

Po,m(<x) 2 = G(a)- 1 P , M G'(«)G'(a)- 1 P 0iM G'(a) 



G{a)-'Pl M G{a) 



= P 



0,M 



a). 



(3.7) 



Using this, we have 

Tr S+S-P 0iM (a) 



Tr S+S n P OM {u)Po,M{u) 



< ^Tr Po, M {a)*S+S-S+S~P 0yM {<y) ■ Tr P 0iM (a)*Po,M(a) 
< 



o+ rt- 



Tr P ,M(a)*P ,M(a). 



(3.8) 



From the definitions, (3.1) and (3.6), one has 

Tr P ,m(«)*P ,m(«) = Tr G(a)P , M G'(a)- 1 G'(«)- 1 P , M G'(a) 

= Tr G(a)- 2 Po, M G(a) 2 P ,M 
< q M \\PoM 2a )\\ ■ 

Combining this, (3.3), (3.5) and (3.8), we have 



(3.9) 



m n 'Q,M 



< 



P 0i m(2o;)|| exp —a (l - R- K ) {dist(m, n)} 1 -^ 2 ] . (3.10) 



Therefore it is sufficient to evaluate the norm ||P ,Af(2o;)|| with a suitable choice of the 
parameter a. 

Let us introduce the contour integral representation of the projection P 0jM as 



Po. 



M 



1 r dz 
2m Jr z — Pa 



(3.11) 



where the closed path T is taken to encircle all of the eigenvalues in the sector of the 
corresponding ground state, and Pm is the projection onto the eigenspace H,k,m of S\ 
with the eigenvalue M. From the definition (3.6), one has 

1 t dz 

P , M (2a) = G(2a)- 1 — / — P M G(2a) 

Zm Jr z — H\ 

dz 



dz 



2m Jr z- H' A 



M 



(3.12) 



with the non-hermitian matrix H' A = G{2a)~ 1 H^G{2a). In order to evaluate the resolvent 
(z—H' K )~ l in the right-hand side, we begin with getting the explicit form of the transformed 
Hamiltonian H' k . Note that the Hamiltonian °f (2-5) is written as 



{i,j}eA b 



(3.13) 
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Using the definition (3.3) of and the relations (3.4), one has 
G(2a)- l Hl Y G(2a) = 



{i,j}eA b 



Hf Y + K A + %L t 



with the two hermitian matrices, 



K A = £ J™{cosh[2a(9 i -9 j )]-l}(S?S- + S-S+ 

{i,j}€M,R(m) 



and 



L A = -i E J ?J sinh[2a(^ - 9 3 )} (Sf Sj - S~S] 

where the set A ltR (m) of the bonds is given by 

^i,fi( m ) — {{hj} e j 1 < dist(i,m) < i? — 1, dist (j, to) > dist(i,m)}. 
Thus we have 

H' A = H A + K A + iL A . 
The norms of the operators, K A and L A , can be estimated as follows: 

Lemma 3.1 The norm of the operator K A satisfies 



(3.14) 
(3.15) 
(3.16) 

(3.17) 
(3.18) 



I^aII < J, 



XY 



(SiSj + S.S+)\\cl{ C osh2a-l) 2 ^ +D 1 



+ D-2 



(3.19) 



Remark: The bound implies that one can make the contribution from K A small in the 
resolvent (z — H'^ 1 by choosing a small a. 
Proof: From (3.15), one has 

II^aII < Jl Y x \\{stSj + STS})\ x E {cosh[2a(0i - 0,-)] - 1}. (3.20) 

{i,j}eA ltR (m) 

The sum in the right-hand side is rewritten as 

E {oosh[2a(fl i -fl i )]-l}= E E E {cosh[2a(^.)]-l}. (3.21) 

{i,j'}6Ai ji{ (m) r— 1 i:dist(i,m)=r i:{j,j}GA 6 

dist(j,m)=r+l 

From the definitions (3.2) and (3.3), one has 

\6i -6j\< kET^ 1 - ' 2 ^- 1 < 1 (3.22) 



for i,j satisfying dist (i, to) — r and dist (j, to) = r + 1 for r = 1, 2, . . . , R — 1. From the 
second inequality, 

rnsU'Wfl. - OA] — 1 rnsh9rv-1 

(3.23) 



cosh[2a(0i - 0,)] - 1 < cosh 2a - 1 



4a 2 (9i -9j) 2 



4a 2 
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Combining these inequalities, one obtains 

cosh[2a(0i - 6j)\ - 1 < (cosh 2a - 1)(^ - 6j) 2 < (cosh 2a - l)iT 2K+2 " D r 2 ^ 2 . (3.24) 
Substituting this into the right-hand side of (3.21), we have 

£ {coshpo^ - 9j)] - 1} 

{*j}GAi i _ R (m) 

< (cosh2a-l)X; X! E i?- 2K+2 - D r 2K - 2 

r=l i:dist(i,m)=r Aj, 
dist(j,m)=r+l 

< C7 (cosh2«-1) J2 ]T R- 2K+2 ~ D r 2K - 2 

r=l i:dist(i,m)=r 

< C 2 (cOsh2a - 1) J2 r D-l R -2n + 2-D r 2 K -2 

r=l 

2k 4- D — 1 

< C7 2 (cosh2 Q -l) 2/t + D _ 2 , (3.25) 

where we have used the assumption (2.2) on the (fractal) dimension D and the definition 
(3.2) of the parameter k. Substituting this into the right-hand side of (3.20) gives the 
desired bound (3.19). ■ 

In a similar way, we can obtain the following bound for \\L\\\: 

Lemma 3.2 The norm of the operator La satisfies 

R D / 2 \ 



(S+Sj - SrSf)\ C 2 \ sinh2a| ^1 + 



T II < T XY 

A II ' 'J mav I ' i I ill' ■ — ' i y-. 

K + D — 1 



(3.26) 



Proof: In the same way as in the proof of Lemma 3.1, one has 

R-l 



< J XY 

— " max 



{S?Sr-SrS+)\\xY, E E | sinh[2a(0i - 0,)]!- (3.27) 

r=l i:dist(i,m)=r i:{j,i}GA(, 

dist(j,m,)=r+l 



For $i, 8j in the sum, the following bound is valid: 



sinh[2a(6»i - 0j)]| < | sinh2a| 2g ^ 



|2a(0i-0j)l ~~ l 2 «l 
Therefore we have 

| wrih[2a(6i - 6j)]\ < | sinh2a||0; -6j\<\ smh2a\R- K+1 - D/2 r K - 1 . (3.29) 

Substituting this into the right-hand side of (3.27), we can obtain the desired bound (3.26) 
in the same way as in the proof of Lemma 3.1. ■ 

The upper bound (3.26) increases as the distance R between the two spins increases. In 
fact, this bound is not sufficient for estimating the resolvent (z — H^)' 1 . We will further 
employ the technique developed in [32, 33]. For this purpose, we need the following 
estimate for the double commutator: 
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Lemma 3.3 The following bound is valid: 

\\[L A ,[H A ,L A }}\\ <C 1 (sinh2 Q; ) 



2 2k + D- 1 



2/t + £>-2' 

where C\ is a positive constant which is independent of the parameter a. 

Proof: Write the Hamiltonian H A in terms of the local Hamiltonian h UtV as 

Ha = E hu,v with h U:V = Ju^(S^S v + S u S^) + J^ V S^S^. 
{u,v}eA b 

Note that 



(3.30) 



(3.31) 



[H A ,L A ] = -t £ E ^ Y sinh[2a(^-^)]K„(^7- 5 r^ 

{«,t)}GA[, {i,j}GAi ifl (?Tt) 

= -i E j5 Y sinh[2a(^-^-)] E [^(SfST-SrSf) 

{i,j}eA 1}R (m) {u,v}n{i,j}^9 

-i E ^sinh[2a(^-^)] E M ^ (3-32) 

{«J}6^i,flH teA s :dist(t,{i,j'})=0,l 



where M ij;t is a matrix with the support {i,j, t} C A s . Using this, the double commutator 
is written as 

[L A ,[H A ,L A }\ = - E E J fc X 7jJ Y sinh[2a(^-^)]sinh[2a(^-^)] 

{k,£}eA ltR (m) {i,j}eA hR (m) 

E [(^"-S^.MyJ. (3.33) 

teA a :dist(t,{i,j})=0,l 

In the same way as in the proof of Lemma 3.1, we have 



\\[L A ,[H A ,L A ]]\\ < (J*lf (sinh2a) 2 J R- 



2n+2-D 



„f\K— l„/t— 1 



_R-1 

x E E E E E E (O-r 

r'=l fc:dist(fc,m)=r' {fc,£}eA 6 r=l i:dist(i,m)=r {jj'leA;, 

^:dist(£,m)=r' + l j:dist(j,m)=r+l 

7- c + 



E 

teA s :dist(t,{i,j})=0,l 



(S k S e - S k S e ), Mij. t 



(3.34) 



We decompose the sum in the right-hand side into two parts, I± with r > r' and I 2 with 

r' > r, as 



„/\k— i„«— i 



J i = E E E E E E Wr r 

r'=l k:dist(k,m)=r' {k,£}eA b r>r' i:dist(i,m)=r {i,j}eA b 

£:dist(£,m)=r'+l j:dist(j,m)=r+l 

E 

teA s :dist(t,{iJ})=0,l 



(3.35) 
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and 



= E E E E E E (r'rv- 1 

r=l i:dist(i,m)=r {ij'jeA;, r'>r k:dist(k,m)=r' {k,£}eA b 

j:dist(j,m)=r+l £:dist(£,m)=r' +1 

E 

t6A s :dist(t,{j,j})=0,l 



(3.36) 



First let us estimate J2. Since r' > r in the sum, one has (r') K 1 < r K 1 . Using this 
inequality, I 2 is evaluated as 



„2k-2 



R-l 

h < E E E 

r=l i:dist(i,m)=r {«,j}eA(, 

j:dist(ji,m,)=r+l 

x E E E E 

r'>r fc:dist(fc,m)=r' {fe,^}GA 6 *6A s :dist(t,{i,j'})=0,l 
e-.dist(e,m)=r'+l 

< E E E - 2K " 2 

r=l i:dist(i,m)=r {«,j}sA(, 

j:dist(j,m)=r+l 

E E |[(^--^/),M M;i 

teA s :dist(t,{ij})=o,i {M}n{i,j,t}^0 

< 12C 4 max{||[(^- - 5^),MyJ||} E^V 



)-l J2.K-2 



r=l 



< 12C 4 max{||[(^- - S^Sf), M hJ , t ] ||} 



1 + 



2k+D-2 



2k + L> -2 



(3.37) 



Similarly we have the same upper bound for I 1 . Combining these bounds, (3.34), (3.35) 
and (3.36), the desired result (3.30) is obtained. ■ 

Now let us estimate P ,m(2q;) of (3.12). The contour integral in the right-hand side is 
written 



H" / 



idy f E 
Je 



dx 



yo E + + iy- H A 
y° idy 
E_+iy-H' A 



M 



e+ x + iy -H' A 

Pm + / : tttPm- 

Je- x-iy - H' A 



(3.38) 



Here we choose the three real numbers, yo, E + and E_, as follows: Relying on Lemma 3.2, 
we take 

y = C 2 R D/2 (3.39) 



satisfying 



C 2 R D/2 - \\L A \\ >C 3 >0 



(3.40) 
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with some positive constants C 2 and C 3 . From the assumption on the spectrum of the 
Hamiltonian if A , we can take E + , E_ so that the distance between the spectrum and 
{E + ,EJ\ is greater than or equal to AE/2. From (3.38) and (3.39), we have 

||iW2a)|| < — R D/2 \ sup \\R'(E + + iy)P M \\ + sup \\R'(E_ + iy)P M \\ \ 

+ E+ ~ E ~ I sup \\&( x + iy )\\+ sup \\R'(x-iy )\\\, (3.41) 

*K [xe[E-,E+] xe[E-,E+] J 

where we have written R'(z) = (H' A — z)^ 1 . 

In order to estimate the norm of the resolvent R'(z) with z = E + + iy, we employ the 
technique developed in [32, 33] for the following matrix element ($ + , L A <£>_): 

Lemma 3.4 Let $_ = P ,m$ and $+ = (1 - P ,m)$ for a vector $ e Ha,m- Then we 
have 

|($ + ,L A $_)|</(a,A)||$ + ||||$_||, (3.42) 

where 

/(«, A) = ^-^-\\[L A ,[H A ,L A ]}\\+2^\\L A r (3.43) 
with AS = max /ij/1 /{|£'o )/1 - E ^\}. 
Proof: Note that 

|($ + ,L A $_>| 2 = ($„,L A $ + >($ + ,L A <l>_> 

< ||$ + || 2 ($_,L a (1-P ,m)£a$->- (3.44) 

The matrix element in the right-hand side is evaluated as 

($_,L A (l-P 0iM )L A $_) < -^($_,L A (H A -E )(1-P , M )L A $. 

>_,L a (# a -£ )£a$- 



AE 

L A (H A - E )P , m La$- 



AE 

< ^($_,L A (iJ A -P )L A <l._) + ^||L A || 2 ||<l._|| 2 , 

(3.45) 

where we have written Eq = X)«=i Eq^/qm- Further the first term in the right-hand side 
in the last line can be written 



$_,L a (# a -Po)£a<£-) = ($-,L A [# A ,L A ]$-) + ($-,L A (# A -£ )$-). 

(3.46) 
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Therefore the matrix element can be evaluated as 

;$_,L A (ij A -Eo)^A$- 

< ^\(^_,L A [H A ,L A ]^) + (L A [H A ,L A ]^_,^_)\ + A£\\L A \\^4 2 
1 



< 



-\(^,[L A [H A ,L A ]]^.)\+A£\\L A \\ 2 \\^. 
l -\\[L A [H A ,L A ]]\\ + AS\\L A \\ 2 ) ||$_f. 



(3.47) 



Combining this, (3.44) and (3.45) gives the bound (3.42) with (3.43). 



By using this lemma and the improved Combes-Thomas method [31], we obtain the 
following lemma: 

Lemma 3.5 Let z = E + + iy with y G R. For a sufficiently large volume \A S \, there exist 
positive constant, a and C 4 , such that 



'H' - z)- l P, 



M 



< C4 for any a < ckq. 



(3.48) 



Both of the constants, a and C 4 , are independent of the volume \A S \ and of the distance 
R between the two spin operators, S+, S~ . 



Proof: Using the Schwarz inequality, one has 

||ct>|| \\(H' A - z)$\\ > Re (($+ - $_), (H' A - z)(Q + + $_)> 



(3.49) 



for any vector $ G H AtM , where $ + = (1 - P ,m)$ and $_ = P ,m$- 

We recall the expression, H' A = H A + K A + iL A . For the hermitian part of H' x — z, one 

has 



Re (($+ - $_), (H A + K A — £+)($+ + $_)) 
> (E 1 - \\K A \\ - E + )\\$ + \\ 2 + (E + - \\K A \\ - max{£ ,,})||<l>-|| 2 

+ Re(($ + ,ir A $_)-(<l>_,K A <l> + )), 



(3.50) 



where E\ is the energy of the first excited state, and Eo ifl are the eigenvalues of the ground 
state. Since the last term in the right-hand side is equal to zero, one has 

Re (($+ - $_), (H A + K A - £+)($+ + $_)) > Qa£ - \\K A \\^ ||$|| 2 , (3.51) 

where we have used the fact that the distance between the spectrum of H a> m and the 
energy E + is greater than or equal to AE/2. 
For the rest of H' k — z, one has 

Re <($+-$_),(iL A -*?/)($+ + $-)) = -Im(($ + -$_),(L A -</)($+ + $_)) 

= -rm(($ + ,L A $_)-($_,L A $ + )) 
= -2 Im($ + ,L A $_) . (3.52) 
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From (3.42), (3.49), (3.51) and (3.52), one obtains 

||$||||(# A -z)$|| > (\ AE ~ II^a||) ||<f|| 2 -2/(a,A)||$ + ||||<l>_ 

-AE- \\K A \\-f(a, A) 
-AE-\\K A \\-f(a,A) 



> 



$|| 2 + /(a< A)(||$ + ||-|I<M) 2 
$|| 2 . (3.53) 



From the assumption (2.3) on the quasidegeneracy of the ground state, Lemma 3.3 and the 
expression (3.43) of f(a, A), one can find that f(a, A) becomes small for a small parameter 
a, and for a sufficiently large volume |A S | compared to the distance R between the two 
spin operators. Combining this observation with Lemma 3.1, we have that there exist 
positive constants, a and C4, such that 



-AE-\\K A \\-f(a,A)>C, 



(3.54) 



for any a < ao, and for a sufficiently large volume |A S | compared to the distance R between 
the two spin operators. Substituting this inequality into the right-hand side of the above 
bound (3.53), we obtain 

\\(H' A -z)$\\ > C 4 \m. (3.55) 
Choosing <3> = (H' A — z)~ l P M ^> with a vector ^, we have 



\P M M\ > C 4 



(3.56) 



Similarly one can obtain the following two lemmas: 
Lemma 3.6 Let z = E_ + iy with y G R. There exists a positive constant C 5 such that 



[H' A - z) 1 P M < C 5 for any a < a , 
where a is the same as in the preceding lemma. 
Proof: Using the Schwarz inequality, one has 

||<t>||||(//;-z)<t>|| > Re^,(H' A -z)^) 

> (nun{^} - H^aII - £_) 

> QaE-||^ a ||)||<1.|| 2 



l$l 



(3.57) 



(3.58) 



for any vector $ G 7Ya,m- Here we have used min^l-Eo^} — E_ > AE/2. Therefore, in 
the same way as in the proof of the preceding lemma, one can prove the statement of the 
lemma. ■ 
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Lemma 3.7 Let z = x ±iy with x G R. Then 



(3.59) 



Proof: Using the Schwarz inequality and the definition (3.39) of yo with the condition 
(3.40), one has 



|ct>||||(#;-z)<t>|| > \Im{<t>,(H' A -z)<*>)\ 
= |($,(L A =F2/o)$>| 

> |(yo-||^A||)|||$|| 2 >C 3 ||<l>|| 2 



(3.60) 



for any vector Taking $ = (if A — z) l ^ with a vector the desired bound can be 
obtained. ■ 

Combining the bound (3.41) with the three Lemmas 3.5, 3.6 and 3.7, we have 

||Po,M(2a)|| < CR D/2 + C for any a < a (3.61) 

with the positive constants C and C. Substituting this into (3.10) and choosing a = ao, 
we obtain the bound (2.10) for the spin-spin correlation. 
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